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Abstract 



> 

' We derive the radiation pressure force on a non-relativistic moving plate in 1+1 di- 

OO ■ mensions. We assume that a massless scalar field satisfies either Dirichlet or Neumann 

, boundary conditions (BC) at the instantaneous position of the plate. We show that when 

I the state of the field is invariant under time translations, the results derived for Dirichlet 

and Neumann BC are equal. We discuss the force for a thermal field state as an example 
for this case. On the other hand, a coherent state introduces a phase reference, and the 
CL|' two types of BC lead to different results. 



1 Introduction 

An intriguing feature of the static Casimir effect is the dependence of the force on the type 
of boundary conditions imposed on the field. For two parallel infinitely permeable plates the 
force turns out to be attractive, and in fact identical to the original Casimir result for a pair 
of perfectly conducting plates. On the other hand, the Casimir force between a perfectly 
conducting plate and a permeable one is, surprisingly, repulsive j2]. Due to this peculiarity, 
permeable plates have been considered recently in the literature in the context of Casimir effect 
as well as Cavity QED "H]. 

An analogous situation takes place for a scalar field in 1+1 dimensions: taking Neumann 
conditions at the boundaries results in the same attractive force obtained with Dirichlet BC, 
whereas in the mixed case the force is repulsive jl]. For reviews on the Casimir effect see [3] 
and references therein. 

In the case of moving boundaries |6J, the force exerted by vacuum fluctuations usually 
contains a dissipative component |f| . The amount of mechanical energy dissipated is converted 



1 



into pairs of real particles . These radiation reaction forces on moving bodies may appear even 
in the case of only one moving wall [Oj^fl]- Connections between the dynamical Casimir effect 
and several interesting phenomena such as the Unruh-Davies effect ^21; black hole physics [T3] . 
sonoluminescence mass corrections J3j and quantum decoherence ^Hl have been discussed 
in the literature 

In this paper, we extend the discussion of the role played by different BC to the dynamical 
effect. We compute explicitly the radiation pressure force for both Dirichlet and Neumann BC, 
and discuss the class of field states for which the force is the same for these two types of BC. We 
consider a single moving boundary in the nonrelativistic approximation, and a massless scalar 
field in 1+1 dimensions in a general quantum state. This paper is organized as follows: Section 
2 presents the results for Dirichlet BC. In Section 3, we consider Neumann BC and show that 
these two BC lead to the same force when the field state is invariant under time translations. 
As an example of such a state, we discuss the case of a thermal field in Section 4. In Section 
5, we consider a coherent state, which is an example of field state that does not satisfy this 
symmetry. Section 6 is devoted to the conclusions and final remarks. 



2 Dirichlet BC 

In this section, we assume that the field vanishes at the instantaneous position of the plate, when 
measured in the Lorentz frame S' which is co-moving at a given time: (p'ix' ,t')\piate = 0, where 
the prime quantities refer to S'. In terms of the laboratory coordinates, this BC is equivalent 
to P 

0(5g(t),t) = O. (1) 

We solve eq. in the long wavelength approximation, and assume the effect of the motion to 
be a small perturbation. We follow Ford and Vilenkin and define 

(f){x,t) = Mx,t)+6(P{x,t), (2) 

where 0o corresponds to the solution with a static plate at the origin and 6(f) is a small pertur- 
bation which takes into account the effect of the motion of the plate. Hence the unperturbed 
field 00 satisfies the wave equation \I\(f)Q{x,t) = and the BC 0o(O,t) = 0. Its normal mode 
expansion in the half-space x > is^ 



00 (x, t)= I duj\l— sin(c^x) [a^ e-'"* + al e*"*J , (3) 
with 

K,al] = 5{uj-uj'). (4) 

A similar decomposition is written in the half-space a; < 0, with the bosonic operators a^^ and 
their hermitian conjugates aj, replaced by independent operators 6^ and 6^. 

^We are using c = 1 along the text. 
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The field operator 6(f) also satisfies the wave equation \I\5(j) = and is submitted to a BC 
which can be obtained directly from eqs. (0) and (j21) by taking the Taylor expansion around 
X = and neglecting terms of second order in 6q{t) : 

50(O,t) = -5g(t)a,0o(O,t). (5) 

The net force on the moving plate may be written in terms of the suitable component of 
the energy-momentum tensor as follows: 

Fit) = {T'\5q~it),t) - T'\Sq^it),t)) , (6) 

with 

T'\x,t) = \ [{d^fM + (7) 

and 5q^ = 5g ± e (with e 0+). The average (...) is taken over an arbitrary field state. For 
simplicity, however, we assume that the state in the half-space a; > is defined with respect to 
the operators and al exactly as the state for the half-space x < is defined in terms of h^^ 
and 6^. In this symmetrical case, the net force vanishes to zero order of 5q (plate at rest at the 
origin). As an example, if we consider the field to be at thermal equilibrium at temperature T 
at the half-space x > 0, then our assumption implies that the field at the left-hand side of the 
plate is also at thermal equilibrium at the same temperature. 

For Dirichlet BC, the term [dtcj))'^ in does not contribute to the force, and we find 



T}^{8q-^{t),t) = \ 



(9,0o)'(O+,t) + a,0o(O+,t),a,<50(O+,t) \ + 0{5q 



(8) 



where {A, B} represents the anti-commutator of two given operators A and B. 

The Fourier representation allows for an interesting physical interpretation of the dynamical 
Casimir effect. Moreover, it provides a simple ultra-violet regularization of the force (alterna- 
tively, one may employ the point-splitting method in the time domain ^U])- When taking the 
Fourier transform of (jH)), we replace T^^ by the r.-h.-s. of (jHI). As discussed above, the terms 
independent of 5q from each side of the plate cancel in the symmetrical case. On the other 
hand, the linear terms in 5q add to give 

H^) = - j ^({5.$o(0+,a;'),5.5$(0+,a;-^0})- (9) 

The integration variable uj' in (jH)) is the frequency of the unperturbed field $o; we avoid any 
change of variable of integration, so as to conserve its physical interpretation. 

We now solve for the perturbed field 5$ in terms of $o- We take the solution that propagates 
outwards from the plate: 5$(x, u) = e*'^'^'5$(0, cu), and from the Fourier transform of the r.-h.-s. 
of (0) we derive 

/I f 
-^6Q[u-u']d,%[0^,u']. (10) 
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Eq. (fTn|) shows that the scattering of a given field Fourier component uj' generates a new 
frequency (or sideband) at frequency ui' + fi, where VL is the mechanical frequency. We replace 
eq. (Uni) into Q to find 

^dM=^ ! ^{^-uj') [ ^SQ{u-u'-u")aD{uj',u"), (11) 

J Z-K J 271 

where we have defined the correlation function of the unperturbed field operator 

aD{uJ,uj') = {{d,%{0+,u;),d,%{0+,uj')}). (12) 
For future reference, we derive, by taking the Fourier transform of Q, the result 

a^$o(o+,u;) = ^/inhiui e{uj)a^ + e{ -uj)a}_J . (13) 

From ()13|l . we may derive the correlation function aD{^^,uj') and hence the dynamical Casimir 
force for a variety of field states. However, before considering some specific examples, we 
consider in the next section the case of Neumann BC. 

3 Neumann BC 

In this section we assume that space derivative of the field, taken in the instantaneously co- 
moving Lorentz frame, vanishes at the plate's instantaneous position: 

d,'(P'{x',tXiate=0. (14) 

In the nonrelativistic approximation, and using the appropriate Lorentz transformation, this 
BC can be written in terms of quantities in the inertial frame of the laboratory as follows |llj : 



d^ + 6q{t) dtj(j){x,t)\,=sg{t) = 0. (15) 
As in the previous section, we write = 0o + (^0, but now the unperturbed field is given by 

f°° I h 

cj)o{x,t)= rfcuW— cos(cux) [a^e-*"* + aie*"*l , (16) 
Jo \ rruj 

whereas the motion induced correction satisfies, up to first order in 6q{t) and its time-derivatives. 



dMiO^t) = -5q{t)di<Po{M-m) 5t0o(O,t), (17) 
which follows from the Taylor expansion of (|T5|l . 
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The force on the plate is computed from the energy-momentum tensor as discussed in the 
previous section. However, in contrast to the case of Dirichlet BC, here only the term {dt(t)Y 
contribute in the expression of the r.-h.-s. of eq. ((Tj), since (IT^ yields 

{dA?mt).t) = 5q{tf{d,ct>f{5q{t),t) ~ 0{5q{tf) ■ 
Hence, the net force on the moving plate can be written as: 

Fn = ~m<PfiSq-'{t),t) - idt<Pf{Sq-it),t)). (18) 
Keeping only the linear terms in 6q, we obtain: 

FN{t) = -i(|ai0o(O+,t), W(0+,t)| - {0+ ^ 0-}) + 0{5q'). (19) 

As before, the contributions coming from terms involving only the field operator (f)o cancel out, 
and the linear contributions from each side of the plate are equal, so that the net force on the 
plate can be written in the frequency domain as 

:Fn{u;) = I ^(^-a;')^'({'^'o(0+,a;'),5$(0+,^-a;')})- (20) 

It is now convenient to express 6^(0'^, u') in terms of dxS^{0~^ ,uj'), since this last quantity can 
be written with the help of eq. (|17p (after a Fourier transformation) in terms of the unperturbed 
field $0- The outward solution of the wave equation given the value of 8^6^(0, u) is 

piLij\x\ 

6<l>(x,uj) =e(x)dj<!>(0,uj) , (21) 

where e denotes the sign function. Substituting eq. into eq. (jSOI), we obtain: 

:FN{uj) = -t I ^uj'{!^MO+,u'),dJ<l>{0,uj-u')'^) . (22) 
As a last step, we take the Fourier transform of eq. fTTj) and replace the result into (j22I): 

^ / ^{uj-uj')uj'uj"5Q{uJ-u' -uj")aN{oo',Lo") , (23) 
27r J 27r 

where we defined the correlation function of the unperturbed field operator 

ar,{u;,u') = (|$o(0+,cu),<l>o(0+,a;')l). (24) 
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in analogy with (fT^ . The unperturbed field appearing in is computed from (fTH|l : 




$0(0,^) = J-\ \e{iu)a^ + e{-Lu)al] . (25) 



By inspection of eqs. (fT!?|l and ()25|1 . we derive the following relation between the correlation 
functions for Dirichlet and Neumann boundary conditions: 

aoiujjUj') = \uu'\a]\f{u,u'). (26) 

This result allows us to compare the forces for Dirichlet and Neumann boundary conditions. 
When the field state is invariant under time translations, the correlation function (0o(O, t)0o(O, t')) 
is a function oit — t' only, and then in the frequency domain it satisfies a{uj^uj') oc 5{ijJ + uo'). 
Hence we may replace \ujuj'\ = —ujuj' in ()26|). and from and ()23p it follows that Dirichlet 
and Neumann BC provide the same result: 

J^n{u:) = = J'iu:) (27) 

Moreover, since uj' = —u" in (jllj) and (j23|) . J-'{uj) is proportional to 6Q{uj), and may be written 
in terms of a susceptibility function x{^) as follows: 

J^iu)=x{u;)6Q{u). (28) 

In the time domain, (|28|) reads 

F{t)= [ dt'xit-t')5Q{t'), (29) 



where xit) is the inverse Fourier transform of x{^)- Thus, the effect of the plate displacement 
at time t' on the force at time t depends only on t — t', as expected from the assumption of 
time translational symmetry. 

As an example of field state obeying this assumption, we consider in the next section the 
force for a thermal field (temperature T). This example contains the case of the vacuum state 
as the limit T = 0. In Section 5, we compute the force for a coherent state (amplitude a). 
This example also contains the vacuum state as a limiting case (a = 0), but this time we find 
different results for the two types of BC when \a\ > 0. 

4 Thermal state 

We compute in the Appendix the correlation function ctd for a thermal field: 

aD{uJ,uj') = A7fh\uj\{l + 2n{uj))6{uj + uj'), (30) 
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with 

n{uj) = [exp{h\uj\/kBT) - 1]-^ (31) 

representing the average photon number at frequency uj [ks is the Bohzmann constant). As 
expected and discussed in the previous section, aoioj, u') is proportional to 6{lj+u!'), a signature 
of time translational symmetry. Hence the force is the same for Dirichlet and Neumann BC, 
and it is given in terms of a susceptibihty function x{^) according to (j^H|) . From (jllj) and (jHUI) 
we find 



- + n{uj 



(32) 



Thus the susceptibihty is the sum of two contributions: x = Xvac + Xr, where xt is proportional 
to n(|c(j'|) in ()32|1 . whereas Xvac contains the effect of vacuum fluctuations [the term '1/2' in 
(jnH)]. At the zero-temperature limit n = and then x = Xvac- 

Taking a; > (but an analogous argument is valid for < 0) we can extract a finite result 
for Xvac('^) adopting the following regularization prescription: 

Xy!,c[uj) = ih \ -— \uj \[uj + uj ) = ih \im. \ I +/ +/ \ ^—\uj \[uj + uj ). (33) 

J-oo A^oo \J_j^_^ Jq J 2tT 

Since the integrand is odd under reflection around u' = —uj/2, the first and third integrals in 
the rhs of eq. (j^Hjl cancel exactly for any value of A. Therefore, the single relevant contribution 
is the second term, which corresponds to negative frequencies that give rise to sidebands with 
positive frequencies (on the other hand, when considering a negative mechanical frequency uj, 
the sidebands are down shifted, and the contribution comes from the frequency interval [0, —uj])- 
Positive and negative frequencies correspond to annihilation and creation operators, hence their 
mixture is clearly connected to the emission of photons out of the vacuum state. From eqs. ()28|1 
and (jnSl), we find 

Xvac(^) = «^V(6vr), (34) 

or, in the time domain, Fvac = '^'^^l{'^)-> ^ result first obtained in Ref. (TUl, and which 
corresponds to the nonrelativistic limit of the exact expression of Fulling and Davies P] . 
For the contribution of thermal photons, using that n{uj') = n{—uj') we find from ()32|) 

xAco) = r ^ ^[ (35) 

Jq ztt exp[nuj' / kbI ) — 1 

The integral in may be calculated in terms of the Riemann zeta function: 

Xt{uj)=i — UJ. (36) 

Adding up the results of ()34|1 and ()3fi|l we find in the time domain 
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in agreement with Ref. jTHj. By direct comparison of (jH^ and (jHE)), the viscous, thermal 
contribution in ()H7j) is much larger than the vacuum force when the typical Fourier components 
of 6q{t) are much smaller than ksT /h, which is of the order of lO^^s"^ at room temperatures. 

The thermal contribution does not result from the emission of new particles only. It is in 
part an effect of Doppler shifting the frequencies of the incoming thermal photons. The counter- 
propagating photons are up shifted when reflected by the plate by the amount 5uj = 2u6q(t), 
whereas the co-propagating photons are down shifted by the same amount. By momentum 
conservation, in both cases the plate recoils along the direction opposite to its motion. The 
resulting force is proportional to the total power P incident on both sides of the plate: 

i^Doppier = -2P6qit)- (38) 

In the three-dimensional case, P is proportional to (Stefan-Boltzmann law) and so is the 
thermal viscous force JH] , but in the one-dimensional case considered in this paper the thermal 
power is 

Q h ' 

so that the viscous force as given by (jHTj) is twice the Doppler force As discussed in the 
next section, a similar result holds for coherent states with Dirichlet BC. 



5 Coherent state 

The coherent state of amplitude a is defined as an eigenstate of the annihilation operator: 

a^\(y) = a6{uj — u!o)\a), (39) 

where ujq > represents the frequency of the excited mode [20] • The symmetrical correlation 
function for the Dirichlet case is computed from ()13|) and ()39|) : 

aoiuj, uj') = a]^^{uj, to') + a\y{uj, to') + a^(cu, u'), (40) 

with 

a]^'{uj, iu') = ATrh\uj\5{iu + u'), (41) 



(t]j{uj, uj') = 8'Kh\uj\\a\'^ 5{uj — luq) + 6{uj + luq) 5{uj + Uj') 



(42) 



(43) 



a''^ is the contribution of vacuum fluctuations, which are contained in the coherent state \a). 
This term is also present in the thermal case, and corresponds to the zero-temperature limit 
(n = 0) of (jnO))- It originates from the commutator [a(^,a|'^,] as given by (jl)). As expected. 
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<^D = o'Jf^ when a = 0, since \a = 0) is the vacuum state. The normally-ordered correlation 
function is the sum of cx^ and a%. The former originates from terms of the form alja^j', whereas 
the latter originates from terms a^a^^i and their hermitian conjugates. 

We calculate the force for Dirichlet BC from (fTTj) and (pijl - pHjl . and then for Neumann BC 
from and ()26|) . We write separately the contributions from each term in ()40p : 

J'niuj) = x..c{^)SQ{uj) + J'\uj) ± ^"(cu), (44) 

AT 

with Xvac(t^) given by (jHHjl and 

^i(^) = l!^|«|2^5Q(u;), (45) 



a {uj — iUo)SQ{iu — 2iU(j) + a* {u + uJo)6Q{iu + 2iUQ) 



(46) 



As expected, the contribution from ajf'^ leads to the vacuum force already discussed in con- 
nection with the zero-temperature limit of the thermal field [see (j34p ]. The contribution from 
a\) is also of the form JF^I 



X^{uj)6Q{uj) and the same for the two types of BC, because 
a]^{uj,uj') oc 6{uj + uj'). However, the coherent state contains a phase reference, namely 6 in 
|a|e"', which breaks the invariance under time translations. As a consequence, the forces 



a 



for Dirichlet and Neumann BC are different due to the presence of the term JF^^ in (03]). As 
expected, the symmetry is restored if we average over S, i.e. if we replace the coherent state by 
the incoherent statistical mixture defined by the density matrix operator 



P 



2n 



|a|e*^)(|a|e^' 



(47) 



In fact, jF^^(u;) vanishes if we take such an average, as it can be checked in (j46j) or (perhaps 
more easily) from the expression in the time domain: 



I |2 
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cos(2ix;ot — 26)6q{t) + a;osin(2ix;o^ — 26)6q{t) 



(4J 



Since J-'^ioo) does not depend on the phase S, it does not change when the phase average is 
taken, and reads in the time domain 



F\t) 



a\'^6q{t). 



(49) 



It is useful to calculate the power P incident on the plate, in order to compare with the 
Doppler force as given by (j38|) . For this end, we calculate the average of the energy-momentum 
tensor component = {dx(j)i{x,t),dt4>i{x,t)}/2 considering the incident field (j)i{x,t) alone, 
and taken in normal order so as to discard the contribution of vacuum fluctuations: 



.01 ■\ 
1+ •/ 



hujo 



lap + Yleia^e 



(50) 
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for a; > 0, and 




(51) 
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for X < 0. The total incident power is 



P={:T^l:){x = 0-)-{:T^^:)ix = 0^) 



«ni + cos {2uJot-2S)]. 



(52) 
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We find the Doppler force for a coherent state by replacing (j52|l into (j38|) . For Dirichlet BC, 
the sum of the viscous terms in PH|) and ()49j) is twice the Doppler force, as in the case of the 
thermal field for both BC. This also holds for Neumann BC only if we average over the phase 
6. 

For a pure coherent state the force does not vanish when the plate is at rest, due to the 
second term in ()48|) . which is proportional to the instantaneous position of the plate. This is 
a consequence of the spatial dependence of the energy-momentum tensor of the incident field 
[see (|5U|) and (|5ip]. On the other hand, for the thermal state and the incoherent mixture given 
by ()47p . translational symmetry forbids the existence of a net average force on a single static 
mirror. 

6 Conclusion 

Dirichlet and Neumann BC yield the same force on a moving mirror (in the nonrelativistic 
approximation) when the field state is symmetrical under time translations. This is the case 
for vacuum and thermal states, but not for a coherent state. The coherent state introduces a 
phase reference which breaks this symmetry. When averaging over the phase of the coherent 
state, we recover identical results for the two BC. 

For both thermal and coherent states, the force may be split into two contributions: the first 
accounting for vacuum fluctuations, the second representing the contribution of the normally- 
ordered correlation function. The latter vanishes in the limits of zero temperature (for the 
thermal state) and of zero amplitude (for the coherent state). Thus, the force for the vacuum 
state may be recovered from these two examples as a limiting case. 

In both thermal and coherent cases the force contains a dissipative component proportional 
to the velocity of the mirror (viscous force). It is tempting to interpret this effect as a con- 
sequence of Doppler shifting the frequencies of the incident photons. However, the change of 
amplitude and phase by reflection is also determinant. In particular, the phase acquired by 
reflection underlies the difference between the results for Dirichlet and Neumann BC. No ex- 
plicit connection between the viscous force and the Doppler effect holds for the coherent state 
with Neumann BC, whereas in the other examples discussed in this paper the viscous force was 
found to be twice the Doppler force. 

C.F. and P.A.M.N. thank CNPq for partial financial support. P.A.M.N. also acknowledges 
support by PRONEX and the Brazilian Institute for Quantum Information. 
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APPENDIX: Correlation function for a thermal state 

In this Appendix we compute the correlation function (70(00,00') for a thermal field, used in 
Section 3 for deriving the force. From eq. (jl3|) we find 

aoioo^oo') =C{oo,oo') +C{oo',oo) (53) 

with 

C{oo,oo') =A7ihy/oo\oo'\ e{oo)e{-oo'){a^al^,) + AttH ^y\oo\oo' e{-oo)e{oo'){al^a^>). (54) 
Using the thermal correlations 

{au^alj = (fi{oo) + 1^5{oo + oo') 

{a-ujauj') = n{ij') 6{lj + uj'), 
where n{oo) is the average thermal photon number defined in (jHT|) . we find 



aniuj, 00') = Anh 



ooe{oo) 1 + 2n{oo) + oo'e{oo') 1 + 2n{oo') 



6{oo + oo'). (55) 



Finally, we use that n{—oo) = n{oo) and uj{p{uj) — 9{—oo)) = \ou\ to derive pUj) from 
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